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ABSTRACT 

We present an analysis of the 3-79 keV NuSTAR spectrum of the low-luminosity active galac¬ 
tic nucleus NGC 7213. In agreement with past observations, we hnd a lower limit to the 
high-energy cut-off of > 140 keV, no evidence for a Compton-reflected continuum, and 
the presence of an iron Ka complex, possibly produced in the broad-line region. From the 
application of the MYTORUS model, we find that the line-emitting material is consistent with 
the absence of a significant Compton reflection if arising from a Compton-thin torus of gas 
with a column density of S.Ol^'g x 10^^ cm“^. We report variability of the equivalent width 
of the iron lines on the time-scale of years using archival observations from XMM-Newton, 
Chandra and Suzaku. This analysis suggests a possible contribution from dusty gas. A fit with 
a Comptonization model indicates the presence of a hot corona with a temperature kT^ > 40 
keV and an optical depth r < 1, assuming a spherical geometry. 

Key words: galaxies: active - galaxies: Seyfert - X-rays: galaxies - galaxies: individual: 
NGC 7213 


1 INTRODUCTION 

The central engine of low-luminosity active galactic nuclei 
(LLAGNs) is thought to be powered by accretion of surrounding 
matter on to a supermassive black hole, similar to more powerful 
AGNs, like Seyfert galaxies and quasars (see |Ho| 12008] for a 
review). The X-ray spectrum of AGNs is generally dominated by a 
primary power law component, which is thought to be produced by 
Comptonization of optical/UV photons emitted by the underlying 
accretion disc in a hot plasma, the so-called corona (see, e.g., 
|Haardt & Maraschi|1991[[Haardt et al.|1994||1997^ . A signature of 

this process is the presence of a high-energy cut-off in the X-ray 
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emission, which has been observed in a number of sources (see. 


e.g., Perola et al. 

2002} 

Malizia et al.|2014| 

IBrenneman et al.|2014l 

Marinucci et al. 

|20T4 

Ballantyne et al.| 

2014| Balokovic et al.| 

20151 lUrsini et a 

□2015 



The distinctive characteristic of LLAGNs is their intrinsic 
faintness (Lboi < 10^'* erg s“^). Moreover, the mass accretion rate 
of LLAGNs is generally small; in terms of the Eddington ratio, 
most of them have L/L^a < 10“^ while luminous AGNs have 
LjL em ~ 0.01 — 1 (see, e.g., Panessa et al. 2006 Kollmeier et al. 


|l. Whether LLAGNS are simply a scaled-down version of 
classical AGNs is a matter of debate. A standard geometrically 
thin, optically thick accretion disc ( |Shakura & Sunyaev||1973] l 
may power LLAGNs jMaoz|2007^, as is commonly assumed for 
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luminous AGNs. However, radiative-inefficient accretion flows 
(RIAFs; see, e.g., |Narayan & Yi||T994^ have been proposed to 
explain some observational properties, in particular their lack of 
a UV bump (see, e.g., |Ho|2009[ |Yu et al.||2011^ . Furthermore, in 
luminous AGNs the hard X-ray photon index and the Eddington 
ratio are positively correlated (see, e.g., [Sobolewska & Papadakis 


Cao 


2009 1 , while an anticorrelation is found in LLAGNs ( |Gu &| 


2009^. This result is consistent with the X-ray emission 


of LLAGNs originating from Comptonization in RIAFs, and it 
suggests a similarity between LLAGNs and black-hole X-ray 
binaries (BHXRBs) in the low/hard state, where the accretion rate 
is low ( |Wu & Gu]|2008) . The low accretion rate might not only 
affect the structure of the inner accretion flow, but also that of 
the putative obscuring torus. The torus is predicted to disappear 
in low-luminosity sources by models depicting it as a clumpy 
wind, arising from the outer accretion disc (see, e.g., |Elitzur 
|Shlosman||2006| and references therein), or as a collection of 
many self-gravitating, dusty molecular clouds (see, e.g.,|Hbnig &| 
|Beckert|2007) . 

NGC 7213 is a nearby {z = 0.005839, as given in the NASA/IPAC 
Extragalacti c Database; distance 25.80 Mpc, IPereira-Santaella 

^1 Oni nl T T A/^XT T _ 1 'T vy 


et al. 20 lO) LLAGN with Lboi = 1.7 x 10*^^ erg s~^ l Em- 
|manoulopoulos et al.| [2012^ that hosts a supermassive black 
hole of ~ 10“ solar masses (estimated from the stellar velocity 
dispersion, see |Woo & Urry||200^ , yielding an Eddington ratio 
of 1.4 X 10“®. It has been historically classified as a Seyfert 
1 because its optical spectrum shows broad emission lines, i.e. 
with a full width at half maximum (FWHM) of a few thousand 
km/s jPhillips|[T?7^ . However, it has also been classified as a 
low-ionization nuclear emission region galaxy (LINER) because 
of the low excitation observed in the narrow-line spectrum jFilip-| 
Ipenko & Hal p em|1984| l. Low-ionization lines were also detected 
in the soft X-ray band ( [Starling et al.|[2005| l. jWu et ah] jl983^ 
measured a UV flux higher than the extrapolated optical flux, 
thus indicating a possible UV bump, but it is still weaker than in 
most Seyfert galaxies. More recently, [Starling et al.| ( [2005| ( found 
no evidence for an optical/UV bump using XMM-NewtonlOM data. 


The X-ray spectrum of NGC 7213 shows peculiarities as 
well. The Compton reflection component is found to be weak or 
absent with XMM-Newton, BeppoSAX and Suzaku observations 
(see, e.g., [Bianchi et al.[[200^ [Lobban et al.[[20T^ , in contrast 
to what is commonly observed in Seyfert 1 galaxies. This may 
suggest that the accretion disc is truncated in the inner region, 
perhaps replaced by a Compton-thin RIAF ( [Lobban et al.[[20l0l (. 
An iron line complex is clearly detected between 6.4 and 7 
keV, consisting of three narrow Ka emission lines respectively 
from neutral Fe, Fe XXV and Fe XXVI (see, e.g., the analysis 
of Chandra data by [Bianchi et al.[|2008) . Given the lack of an 
observed Compton reflection hump, such lines cannot originate 
from Compton-thick material, like the accretion disc or a parsec- 
scale torus. [Bianchi et al^j2008) found the FWHM of the neutral 
Fe Kq line to be consistent with that of the Hq line (~ 2600 
km/s), thus suggesting a common origin in the broad-line region 
(BLR). Another possibility is reprocessed emission from dusty 
gas, which can produce strong, neutral Fe Ka emission with 
a weak associated continuum ([Gohil & Ballantyne||2015|. [Em-[ 
[manoulopoulos et al.[ ( [2012l l found an anticorrelation between the 
X-ray photon index and the X-ray luminosity that, together with 
the low accretion rate of the source, indicates a spectral behaviour 
similar to that of BHXRBs in the hard state. Finally, [Bianchi[ 
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Figure 1. The Vm 57]4/?/FPMA and FPMB 3-79 keV light curve. Bins of 
1000 s are used. 


lower limit of 3501 




In this paper, we report on a NuSTAR observation of NGC 7213 
performed in October 2014. The primary focus of this work 
is modelling the broad-band X-ray spectrum of the source and 
constraining the origin of the X-ray emission. In Sect. 1^ we 
describe the observations and data reduction. In Sect. |3 we 
present the analysis of the 3-79 keV spectrum, fitted with both a 
phenomenological and a realistic Comptonization model, and we 
study the time evolution of the iron line complex by analysing 
archival XMM-Newton, Chandra and Suzaku data. In Sect. we 
discuss the results and summarize our conclusions. 


2 OBSERVATIONS & DATA REDUCTION 

NuSTAR jHarrison et al.[[2013^ observed NGC 7213 starting 
on 2014 October 5, with a net exposure of 109 ks (Obs. Id. 
60001031002). The NuSTAR data were reduced with the standard 
pipeline (NUPIPELINE) in the NuSTAR Data Analysis Software 
(NUSTARDAS, vl.3.1; part of the HEASOFT distribution as of 
version 6.14), using calibration files from NuSTAR CALDB 
v20141107. Spectra and light curves were extracted from the 
cleaned event files using the standard tool NUPRODUCTS for each 
of the two hard X-ray telescopes aboard NuSTAR, which have 
corresponding focal plane modules A and B (FPMA and FPMB). 
The spectra from FPMA and FPMB are analysed jointly, but are 
not combined, allowing for a free cross-calibration constant. The 
source data were extracted from circular regions (radius 75 arcsec), 
and background was extracted from a blank area close to the 
source. Finally, the spectra were binned to have a signal-to-noise 
ratio (SNR) greater than 3 in each spectral channel, and not to 
over-sample the instrumental resolution by a factor greater than 2.5. 

In Fig. [T] we show the AuSTAR/FPMA and FPMB light curves 
obtained in the 3-79 keV energy range. The total variation in the 
light curve is of the order of a few per cent. Therefore, in the 
following we analyse the time-averaged spectrum of the source. 
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3 SPECTRAL ANALYSIS 


Spectral analy sis and model fitting was carried out with the XSPEC 
12.8 package l Arnaud|l996| , using the minimisation technique 
throughout. In this work, the errors are quoted at the 90 per cent 
confidence level, if not stated otherwise. 


3.1 The NuSTAR spectrum 

As a first step, we define a phenomenological model, which we 
describe below. We fit the NuSTAR data in the 3-79 keV range, 
allowing for a free cross-calibration constant between the modules 
FPMA and FPMB. The two modules are in good agreement, with 
a cross-calibration factor Ka-b = 0.99 ± 0.01 fixing the constant 
for the FPMA data to unity. 

We modelled the continuum with a cut-off power law, modi¬ 
fied by neutral absorption (PHABS model in XSPEC) from Galactic 
hydrogen with column density Wh = 2.04 x 10^° cm“^ 

|& Lockman||1990y This simple model yields a poor fit (reduced 
X‘‘ = 556/369), with positive residuals between 6 and 7 keV, which 
can be attributed to the known Fe complex between 6.4 and 7 
keV ( [Bianchi et al.|20()3l|2008| l. |Bianchi et al.| ( |2008) found three 
emission lines in the ChandralHEG spectrum, at 6.4, 6.7 and 
6.966 keV. These lines are interpreted as emission from neutral 
Fe, Fe XXV and Fe XXVI respectively. ChandrafHEG data provide 
much higher energy resolution near the Fe complex compared with 
NuSTAR. We thus followed the analysis of [Bianchi et al.H2008^ 
and tested for the presence of three narrow Gaussian lines fixing 
their energies at 6.4, 6.7 and 6.966 keV, leaving the normalizations 
free to vary. The inclusion of a narrow Gaussian line at 6.4 keV 
yields a significantly improved fit (reduced = 417/368). Adding 
a second line at 6.7 keV further improves the fit (reduced x^ = 
388/367), with a probability of chance improvement less than 
4 X 10“^ according to the F-test. The inclusion of a third line 
at 6.966 keV yields a good fit (reduced x^ = 375/367), without 
prominent residuals and a probability of chance improvement less 
than 5 x 10“"^ according to the F-test. We report the results of this 
fit in Table^ and in Fig.|^we show the contour plots of the cut-off 
energy versus photon index. In Fig. we show the data, residuals 
and best-fitting model. To further test the presence of a high-energy 
cut-off, we included the 70-month average Swifl/BAT spectrum 
( [Baumgartner et S^|2013^ , however the results are essentially 
unchanged. 

We tested for the presence of a reflection continuum by in¬ 
cluding the PEXRAV model in XSPEC, which describes neutral 
Compton reflection of infinite column density in a slab geometry 
( [Magdziarz & Zdziarski|1995| ). This model is adequate for reflec¬ 
tion off a standard accretion disc, given the high column densities 
expected for these structures (see, e.g., [Svensson & Zdziarski[ 
|1994|l. We fixed the inclination angle of the reflector to 30 deg, 
appropriate for a type 1 source (e.g., [Nandra et al.[[19^ . The 
photon index and cut-off energy of the incident spectrum were 
tied to those of the primary power law. However, no improvement 
is found, with only an upper-limit to the reflection fraction TZ of 
0.13, at 90 per cent confidence level. In Fig.j^we show the contour 
plots of TZ versus photon index. This is consistent with the iron 
lines originating from Compton-thin material (Ah = 10^^ — 10^® 
cm“^), which does not produce a prominent Compton reflection 
hump jBianchi et al.|2003) . To test this result further, we replaced 
the PEXRAV component and the neutral Fe Kq line with the 


Dickey 



Plioton index 


Figure 2. Primary continuum cut-off energy vs. photon index contours. 
Solid green, red and black lines correspond to 99%, 90% and 68% con¬ 
fidence levels, respectively. 


MYTorus model, which includes Compton reflection and iron 
fluorescent lines from a gas torus with an opening angle of 60 deg 
( [Murphy & Yaqoob|2009[ (. The inclination angle of the torus was 
fixed at 30 deg. The column densities of the scattered and line 
components were linked and free to vary. The normalizations of 
the scattered and line components were tied to the normalization of 
the primary power law, i.e. the standard MYTORUS configuration 
(“coupled” reprocessor; see, e.g., [Yaqoob[|2012) . The assumed 
geometry corresponds to a covering fraction of 0.5. We obtain a 
good fit (reduced x^ = 373/363) and we find a column density of 
g X 10^® cm“^. This result is consistent with the estimate of 
~ 3 X 10^^ cm~^ for the Fe Ka line-emitting material reported in 
[Bianchi et al.[ ( |200^ 

We next replaced the cut-off power law with the Comptonization 
model COMPPS ( [Poutanen & Svensson|1996t in XSPEC. COMPPS 
models the thermal Comptonization emission of a hot plasma 
cooled by soft photons with a disc black-body distribution. We set 
a spherical geometry (parameter GEOM = —4 in COMPPS) for the 
hot plasma and a temperature at the inner disc radius of 10 eV, 
fitting for the electron temperature fcTe and Compton parameter 
y = TT{kTe/meC^). The model yields a good fit (reduced x^ = 
375/366) with best-fitting parameters kTe = 295 I 250 
y — 0.52('(q ®q, which imply an optical depth t = 0.2('(g {. In Fig. 
[^we show the kT^ — y contour plots. 


3.2 Time evolution of the Fe Ka line 

To try to understand the origin of the Fe Ka lines, we compared 
our results with past observations of NGC 7213 by re-analysing 
the archival XMM-Newton (2001, 2009), Chandra (2007) and 
Suzaku (2006) data. We fitted the 3-10 keV data sets with the 
baseline model described in Sect. [^i.e. including a cut-off power 
law and three narrow Gaussian lines. The best-fitting parameters 
are consistent with those reported in [Bianchi et al.[ j2003[ [2008[ l 
and |Lobban et al.[j2010[ l. 

In Fig. 1^ we show the time evolution of the neutral Fe Ka line 
EW and flux, and of the 3-10 keV continuum flux and photon 
index. The continuum flux variations reach up to a factor of 2, 
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Figure 5. COMPPS Compton parameter vs. temperature contours. Solid 
green, red and black lines correspond to 99%, 90% and 68% confidence 
levels, respectively. 


Figure 3. The NuSTAR spectrum and baseline model (see Table^. Upper 
panel: NuSTAR data and folded model. Middle panel: contribution to x^- 
Lower panel: best-fitting model E‘^f{E) with the plot of the additive com¬ 
ponents, namely the cut-off power law and the three Fe Ka lines. 



Figure 4. PEXRAV reflection fraction vs. photon index contours. Solid 
green, red and black lines correspond to 99%, 90% and 68% confidence 
levels, respectively. 



Date (year) 


while the line flux varies at the Icr confidence level only during 
the Chandra/HEG observation of 2007. The line EW, on the other 
hand, exhibits significant variations up to a factor of 1.7. The data 
show a hint of an anticorrelation between the continuum flux and 
the line EW, though we cannot formally establish a significant 
relation. The Spearman’s rank correlation coefficient is —0.4, but 
the test returns a P value of 0.42, implying that the null hypothesis 
is not rejected. Similar results hold for the Fe XXV and Fe XXVI 
lines, albeit with larger uncertainties on their parameters. 


Figure 6. Time evolution of Fe Ko; line parameters and 3-10 keV primary 
flux and photon index. Panel (A): Fe Ka equivalent width in eV units. Panel 
(B): Fe Ko; line flux in units of 10“^ photons cm“^ s“^. Panel (C): contin¬ 
uum 3-10 keV flux in units of 10“^^ erg cm~^ s“^. Panel (D): 3-10 keV 
photon index. Error bars denote the l-c uncertainty. The blue solid lines 
represent the mean value for each parameter, while the blue dashed lines 
con'espond to the standard deviation. 
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Table 1. Best-fitting parameters of the baseline model including a primary 
cut-off power law and three Gaussian lines (line 1 at 6.4 keV, line 2 at 6.7 
keV, line 3 at 6.966 keV). 


r 

Ec (keV) 

^’3-lOkeV ergcm“2 s“l) 

EWi (eV) 

Fi ( 10 “® phcm -2 s-l) 

EW 2 (eV) 

F2 (10-® phcm-2 s-i) 

EW 3 (eV) 

F 3 ( 10 “® phcm -2 s-l) 
X^/d.o.f. 


1.84 ± 0.03 
> 140 
1.3 ±0.1 


98 


+30 

-20 


1 '7-I-0.3 
-*^•'-0.4 
90+24 

2a_ig 

Q r+0.5 

42± 20 
0.6 ±0.3 
375/366 


4 DISCUSSION AND CONCLUSIONS 

We reported results based on a NuSTAR observation of the 
LLAGN/LINER NGC 7213. We derived constraints on the param¬ 
eters describing the high-energy (3-79 keV) spectrum, and studied 
variability over a few year time-scale comparing our results with 
archival observations of the source. 


We have been able to constrain the reflection component, finding 
no evidence for a significant Compton reflection continuum. 
This is consistent with previous results (e.g., |Bianchi et al.|200^ 


[Lobban et al.|201()||Emmanoulopoulos et aH2013| >. We note that 

the constraints on the reflection component are derived assuming 
a static X-ray source ([Magdziarz & Zdziarski|[T995 Murphy & 


Yaqoob|[2009[), but it may actually be outflowing (Beloborodov 


Malzac et al.|200r). An outflowing, energetically dominant 


corona can be generated by a geometrically thin, optically thick 
accretion disc at low accretion rates (?). This scenario is consistent 
with the lack of a significant optical/UV bump from the disc, 
if the power output is dominated by the corona. Furthermore, 
an outflowing X-ray corona could account for the absence of a 
disc reflection component, since the primary emission would be 
beamed away from the disc. This has been suggested especially 
for radio-loud sources having relativistic jets, where the X-ray 


corona might be the base of the jet itself (see, e.g.. 

Walton et al. 

2013 

|Lohfink et al.||2013| |Ballantyne et al.||2014| 

Fabian et al. 


2014 1 . Although NGC 7213 is a radio-intermediate source with 
no clear evidence for a strong jet ( [Blank et al.||2005l [Bell et aT] 
[201 Ijl, its radio emission is found to be weakly correlated with the 
X-ray emission ( [Bell et al.|2011[ l. Therefore, the outflowing corona 
scenario cannot be ruled out in this source. 


We confirm the presence of an iron complex between 6.4 
and 7 keV ( jBianchi et al.|[200^ , consisting of three narrow Ka 
emission lines from neutral Fe, Fe XXV and Fe XXVI. The lack of 
a Compton reflection hump above 10 keV indicates that these lines 
cannot originate from Compton-thick material. Using the MY- 
TORUS self-consistent model, we find that the neutral Fe Ka line 
may be produced by material with a covering fraction of 0.5 and a 
column density of S.O^ti’e x 10^® cm“®. This result is consistent 
with the iron line complex being produced in the BLR |Blanchi| 
jet al.|2008| l. Assuming that the line widths represent a Keplerian 
velocity, it is possible to estimate the distance of the BLR using 
the virial theorem, i.e. Rpc — 4.5 x 10®/“^ Ms (AVkm/s)~^, 


where i?pc is the distance in units of parsecs, / is a coefficient 
that accoun ts for the unknown geometry and orientation ([Peterson] 


et al. 2004 1 , Ms is the black hole mass in units of 10 ° solar masses 
and A Vkm/s is the line width in units of km s“®. [Bianchi et al.[ 
(2008|l report a FWHM of 24001^5]° km s"® for the neutral 
Fe Ka line using the ChandrafHEG data, and 2640lgJ° km s“® 
for the Ha line, using ESO-NTT optical data. Setting / = 4.3 
( [Grier et al.[[2013) l, the virial theorem yields a BLR distance of 
~ 0.02 pc. Alternatively, reprocessing may occur in dusty gas, 
since the presence of dust may enhance the neutral Fe Ka EW 
( [Gohil & Ballantyne|2015^ . This effect is due to a suppression of 
the reflected continuum, caused by the reduction of backscattering 
opacity in dusty gas ( [Draine p003l l. This hypothesis is supported 
by the presence of strong silicate emission features in the mid-IR 
spectrum of NGC 7213 ([ Honig et al.[[2010l [Ruschel-Dutra et al.[ 
[ 2014[ l. The dust sublimation radius is estimated to be ~ 0.03 pc 
from the optical luminosity ( [Honig et al.|201()| l. Therefore, if the 
line originates from dusty gas, the line-emitting region should be 
more distant than the BLR, as suggested by the constancy of the 
line flux on a few years time-scale (see below). We also note that 
the lack of any cold reflection from Compton-thick material is 
consistent with the models of clumpy tori, which are predicted to 
disappear in low-luminosity sources jElitzur & Shlosman[[2006] 
[Honig & Beckert|2007^ . 


We analysed the neutral Fe Ka line time evolution over ~ 15 years 
by re-fitting archival XMM-Newlon, Suzaku and Chandra data of 
the source. According to our estimate, the BLR distance is ~ 0.02 
pc, or ~ 20 light days (see also [Kaspi et al.[[2005] [Balmaverde[ 
[& Capettij [2014^ . The observed X-ray variability time-scale 
of NGC 7213 is of several days or weeks ( [Emmanoulo p oulo^ 
[et al.[|2012^ , in agreement with the lack of variability during 
the ~ 100 ks NuSTAR observation. We would thus expect the 
Fe Ka line flux to vary in response to the continuum variations, 
albeit with a time delay due to the light-travel time. However, 
we only detect a marginal variability of the line flux during the 
Chandra observation of 2007. Therefore, a contribution from 
more distant material cannot be ruled out. The Astro-H satellite, 
providing an unprecedented combination of high spectral resolu¬ 
tion and effective area at 6-7 keV, will likely enable us to solve 
the ambiguous origin of the Fe Ka line. For example, the line 
might be the superposition of two components, namely a broad 
component from the BLR and a narrow one produced further away, 
•which Astro-H should be able to disentangle ^Reynolds et al.[20T4] . 

The primary continuum is well fitted by a power law with a 
photon index T = 1.84 ± 0.03, and an extrapolated 2-10 keV 
flux F 2-10 = (1.6 ± 0.1) X 10“®® erg cm“® s“®. The 3-10 keV 
continuum flux and photon index that we found for the archival ob¬ 
servations do not show a clear trend (see Fig.|^. The extrapolated 
2-10 keV luminosity is I/ 2-10 = (1.2 ±0.1) x lO'®® erg s“®. Using 
the 2-10 keV bolometric correction of [Marconi et al.[ ( [2004t , we 
estimate the bolometric luminosity to be Lboi = (1.3 ±0.1) x 10®® 
erg s“®. For a black hole mass of 10® Mq, the Eddington luminos¬ 
ity is 1.2 X 10®® erg s“®. These estimates yield an accretion rate of 
~ 0.1 per cent of the Eddington limit, in rough agreement with the 
value of 0.14 per cent by [Emmanoulopoulos et aL|]2012[ l based 
on the broad-band spectral energy distribution. An Eddington ratio 
~ 10“® is lower than what is typically found for luminous AGNs, 
namely ~ 0.01 — 1, and it lies within the theoretically predicted 
RIAF regime ( ^Ho[2009| l. 
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We can only place a lower limit on the presence of a high- 
energy cut-off, Ec > 140 keV, consistent with the lower limit of 
350 keV found by |Lohban et al.| ( |2010| > using Suzaku+SwiftlBAH 
data, and marginally consistent with the value of keV found 

by |Bianchi et al.| ( |2003| l using XMM-Newton+BeppoSAX data. 
Replacing the cut-off power law with a Comptonization model 
and assuming a spherical geometry, we estimate the coronal tem¬ 
perature to be > 40 keV (kTe = 295 I 250 keV) and we measure 
an optical depth r = 0.2to'i. The lack of an upper limit on the 
high-energy cut-off does not conflict with the upper limit on the 
coronal temperature, because a cut-off power law is known to be a 
rough approximation of Comptonization models (see, e.g., |Stern| 
|et al.|[T995) ). First, the high-energy turnover of a Comptonization 
spectrum is much sharper than an exponential cut-off (see, e.g., 
IZdziarski et al.|[2Q03t . Moreover, a Comptonization spectrum is 
actually a superposition of several orders of Compton scattering 
spectra. When the optical depth is small, the different scattering 
orders are separated in energy, thus resulting in a bumpy spectral 
shape (see, e.g., |Poutanen & Svensson|1996| l. But the optical depth 
and temperature are inversely related for a given heating/cooling 
ratio, i.e. the ratio of the power dissipated in the corona to the 
intercepted soft luminosity ( [Stem et al.|1995[[Haardt et al.|1997^ . 
This can produce an upper limit on the coronal temperature even if 
no exponential cut-off is required (see, e.g.,|Petrucci et al.|2Q13^. 


The combination of a weak or absent reflected continuum, a 
weak UV bump, and a low accretion rate suggest that the standard 
optically thick, geometrically thin accretion disc is truncated in 
the inner region of the source ^Starting et al.|2005] [Lobban et al.| 
|2010| >. A possible explanation is that the nucleus accretes via a 
RIAF, with the inner edge of any standard d isc restricted to larg e 
distances, ~ 10^ gravitational radii (see, e.g., Quataert et al. 1999 1 . 
This is a natural suggestion for a LLAGN, since RIAFs are only 
expected in sub-Eddington systems and the low radiative efficiency 
would explain the observed low luminosity. The X-ray emission 
from a RIAF is likely dominated by thermal Comptonization, in 
agreement with observations of X-ray binaries in the hard state 
(see, e.g., |Narayan||2005[ and references therein). The soft seed 
photons can be synchrotron photons produced in the hot accretion 
flow itself (synchrotron self-Compton), or thermal photons from 
the outer thin disc (see, e.g., [Yua n & Zdziarski||2004[ |Nemmeri| 
|et al.|[20T4l ). As we noted above, the X-ray variability time-scale 
of NCC 7213 is of the order of days/weeks jEmmanoulop oulo'^ 
|et al.|[20T^ . Therefore, the Comptonization process is consistent 
with taking place in an extended region with kT^ > 40 keV and 
r < 1, possibly illuminated by the outer thin disc. The relatively 
small optical depth is a further indication that the corona should 
be extended and subtend a large solid angle as seen from the disc, 
in order to scatter a sufficient number of soft photons to produce a 
substantial X-ray continuum. 


Finally, our results can be compared with those on other 
weakly accreting ACNs observed so far by NuSTAR, namely 
NCC 5506 (estimated Eddington ratio as low as 0.7 x 10“^; see 


Matt et al.||2015P and NCC 2110 (estimated Eddington r atio of 


0.25-3.7 xlO ; see Marinucci et al. 2015 1 . Matt et al. j2015| 


^ This value of the Eddington ratio is obtained by assuming the highest 
estimate of the black hole mass in NGC 5506, namely ~ 10® Mq. The 
mass is poorly known, with a lower limit of a few X 10® Mq (see ?), which 
would yield an Eddington ratio above 0.1. 


found a lower limit on the high-energy cut-off of ~ 500 keV for 
NCC 5506, while [Marinucci et al.| ( [^15| l inferred Ea > 210 
keV for NCC 2110. In these objects, then, a weak disc emission 
may be accompanied by a relatively high coronal temperature, 
in agreement with our present work on NCC 7213. The number 
of observed sources with good measurements of the high-energy 
cut-off is too small to establish any statistical correlation with other 
parameters. However, if the X-ray emission is due to a Comptoniz- 
ing corona, a high coronal temperature is expected when the disc 
radiation is weak, because inverse Compton scattering is inefficient 
in cooling the corona. Eurther studies on a greater number of 
sources will be required to confirm this scenario by constraining 
the physical parameters of the disc/corona system. 
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